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Understanding the molecular architecture of the 
flagellum is crucial to elucidate the bending 
mechanism produced by this complex organelle. 
The current known structure of the flagellum has not yet 
been fully correlated with the complex composition and 
localization of flagellar components. Using cryoelec- 
tron tomography and subtomogram averaging while 
distinguishing each one of the nine outer doublet micro- 
tubules, we systematically collected and reconstructed 
the three-dimensional structures in different regions of 
the Chlamydomonas flagellum. We visualized the radial 



and longitudinal differences in the flagellum. One doublet 
showed a distinct structure, whereas the other eight were 
similar but not identical to each other. In the proximal 
region, some dyneins were missing or replaced by minor 
dyneins, and outer-inner arm dynein links were variable 
among different microtubule doublets. These findings shed 
light on the intricate organization of Chlamydomonas 
flagella, provide clues to the mechanism that produces 
asymmetric flagellar beating, and pose a new challenge 
for the functional study of the flagella. 



Introduction 

Flagella and cilia are complex organelles composed of >300 
proteins and are essential for motility and sensory functions 
(Dutcher, 1995; Pazour et al., 2005). Motile cilia and flagella in 
most organisms have the canonical nine microtubule doublets 
(MTDs) plus two central singlets (9 + 2) architecture (Fig. 1, 
A and B). Anchored periodically on MTDs are inner dynein 
arms (IDAs) and outer dynein arms (ODAs), which consist of 
dyneins a-g and a complex containing the a, p, and 7 heavy 
chains, respectively. MTD is connected with the central pair 
through the radial spoke (RS). Dynein arms drive the sliding 
activities between adjacent MTDs, which are then converted 
into the bending motion of the flagellum. The bending activity 
is likely to be regulated by many means that involve the central 
pair and RS (Porter and Sale, 2000), the intermediate chain 
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(IC)/light chain (LC) of inner arm dynein/(dynein II ; Porter and 
Sale, 2000), and the dynein regulatory complex (DRC; Piperno 
et al., 1992), a structure possibly identical with the interdoublet 
(nexin) link (Summers and Gibbons, 1971; Heuser et al., 2009). 
Though sharing these common structures, flagella generate a 
variety of bending patterns. In Chlamydomonas flagella, wave- 
form is asymmetric at low calcium concentrations, whereas 
the waveform turns symmetric at high calcium concentrations 
(Hyams and Borisy, 1978; Bessen et al., 1980). Defects in fla- 
gellar assembly can lead to many diseases, termed ciliopathies, 
such as hydrocephalus, polycystic kidney disease, primary cili- 
ary dyskinesia, and infertility (Fliegauf et al., 2007). Despite a 
long history of research, many questions regarding the molecu- 
lar mechanism of the flagellar function are still to be answered, 
for example, how the bending motion is initiated at the basal 
part of the flagellum or how the flagellum produces both sym- 
metric and asymmetric motions. 
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Figure 1 . Division of the complete Chlamydomonas axoneme. (A) The 

partitions of the axoneme into four regions: PEA, proximal, central, and 
distal. (B) Cartoons and representative 96-nm-thick cross-sections of each 
region from our tomograms, observed from base to tip of the axoneme. The 
markers of different regions are circumferential interdoublet linkers (red) in 
the PEA region, the 1-2 bridge (yellow) in the proximal region, and the 
beaks in MTD1 , 5, and 6 (blue) in the central region. No beak was found 
in the distal region. As a result of the missing wedge effect in the raw 
tomogram, only three circumferential linkers roughly parallel to the optical 
axis are seen. Examination of all our data indicates that the circumferen- 
tial linkers exist between every pair of MTDs in the PEA region, and they 
are absent from other regions. (C) Single tomographic slices showing the 
longitudinal sections of flagella in the central region (top), the PEA region 
(bottom left), and the proximal region (bottom right). The red and yellow 
arrows indicate a short span of the linkers in the PEA region and the 1-2 
bridge, respectively. Bars, 100 nm. 

Structural studies of the axoneme have been contributing 
significantly to tackling those questions. Recent advances in 
cryoelectron tomography (cryo-ET) have allowed researchers 
to obtain the structure of the 96-nm repeating unit of the axo- 
neme at intermediate resolution (^35 A; Nicastro et al., 2006; 
Bui et al, 2008, 2009; Heuser et al., 2009; Movassagh et al., 
2010). However, those reported structures have not yet fully 
reflected the complexity observed by biochemical analysis. 



For instance, according to Chlamydomonas genomics (Merchant 
et al., 2007), more dyneins were predicted in the flagella than 
observed in the 96-nm structural unit within the axoneme. In 
Chlamydomonas, 16 dynein heavy chain (DHC) genes have 
been predicted and classified into five groups: one cytoplasmic 
(DHClb), three outer-arm (ODA a, p, and y), one double- 
headed inner-arm (fa, f/3), and nine single-headed inner-arm 
DHCs (a-e, g, DHC3, DHC4, and DHC11) and an orphan or 
novel DHC (DHC 12) that does not fall clearly into either cyto- 
plasmic or axonemal dynein categories (Wickstead and Gull, 
2007; Yagi et al., 2009). 1 1 axonemal DHCs (a-g, a, (3, and 7) 
can be purified biochemically (Kagami and Kamiya, 1992) and 
localized in the 96-nm unit by EM (Mastronarde et al., 1992; 
Gardner et al., 1994; Yagi et al., 2005; Ishikawa et al., 2007; 
Bui et al., 2008). However, the exact localization of some inner- 
arm dyneins has not yet been achieved because the variety of 
available mutants lacking specific dynein is limited; more spe- 
cifically, the distinction of loci has yet to be made between dy- 
nein a and d and between dynein b and g. In addition, the heavy 
chains of three single-headed dyneins — DHC3, DHC4, and 
DHC1 1 — as well as the orphan DHC 12 (unpublished data) are 
present in small amounts and are therefore termed minor dy- 
neins (MDs). The loci of these dyneins also need to be deter- 
mined. A systematic characterization of the structure of the 
whole axoneme is necessary to pinpoint the location of all the 
axonemal dyneins and, as a result, obtain insight into their inte- 
grated functions. 

So far, all of the 3D reconstructions addressing molecular 
arrangement have assumed identical 96-nm periodic units along 
the entire axoneme and over the nine (or eight) MTDs; this as- 
sumption must be corrected. We previously discovered a radial 
asymmetry in the arrangement of ID As and interdoublet link- 
ages in the central part of the Chlamydomonas axoneme (Bui 
et al., 2009). Seven MTDs have the same dynein architecture, 
which we named the common IDA architecture, whereas two 
MTDs show different structures. The radial asymmetry may 
well underlie the asymmetric planar waveform of Chlamydo- 
monas flagella. With the detection of such a radial asymmetry, 
an obvious next question is whether there are any differences 
along the flagellum. A few studies have dealt with the dynein 
composition specific to the proximal region (Piperno and 
Ramanis, 1991; Yagi et al., 2009). A minor DHC, DHC11, was 
localized in an ^2-um proximal region of flagella, whereas dy- 
nein b was localized outside this region (Yagi et al., 2009). In 
addition, the kinetics of the appearance of DHC3 and DHC4 
during flagellar regeneration also suggested their proximal 
localization (Yagi et al., 2009). Collectively, in the proximal 
region, MDs may replace major dyneins or may be added as extra 
dyneins. To determine this special dynein distribution in more 
detail, structural studies throughout the flagella at molecular 
resolution is necessary. 

In this study, we first tried to determine the arrangement 
of dyneins in the 96-nm periodic unit by analyzing Chlamydo- 
monas wild-type and mutant axonemes by cryo-ET and sub- 
tomogram averaging. We used a novel inner-arm-deficient 
mutant, idalO, which permitted us to distinguish the loci of dy- 
nein a and d as well as those of dynein b and g. Thus, for the 
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first time, the loci for all the major inner-arm dyneins have been 
determined. Then, we systematically acquired tilt series in dif- 
ferent regions of the axoneme while distinguishing each one of 
the nine outer doublet microtubules. This resulted in the finding 
of a structural asymmetry and polarization in dynein localiza- 
tion. These novel features should be important for the studies on 
how the flagellar waveform is determined. 

Results 

Cryo-ET localizes eight: inner-arm DHCs 
in the central region 

Based on the structural indicators in the cross-section (Fig. 1 
[A and B] and Video 1), we divided our tomographic data of 
axonemes into four consecutive regions: proximal end of the 
axoneme (PEA), proximal, central, and distal regions. The 
boundaries between these regions are apparently not sharply 
defined. The PEA region corresponds to an ^300-nm region 
at the very proximal end, which is connected to the transition 
region next to the basal body (Ringo, 1967). It shows no sight 
of the 1-2 bridge, a bridge between MTD1 and MTD2 asso- 
ciated with the proximal region (Fig. 1 B, yellow arrowhead) 
but has circumferential interdoublet linkers (detailed in the next 
section; Fig. 1 C), also referred to as rodlike peripheral links 
in Hoops and Witman (1983) (Fig. 2). The proximal region 
extends for ^2 um and has the 1-2 bridge and the beaks in the 
B-tubules of MTD1, 5, and 6 (Hoops and Witman, 1983; Bui 
et al., 2009). The central region, spanning ^7 um, lacks the 1-2 
bridge. The remaining ^3-um region near the distal end of the 
flagellum, termed the distal region, lacks the aforementioned 
beaks. At the very end of the distal region, there is a short tip 
region with only A-tubule present and in contact with the flagellar 
membrane (Dentler and Rosenbaum, 1977). The structure of the 
tip end was not studied here, as it should be more meaningful to 
study it with intact membrane. 

We reconstructed the 3D structure of each region and 
each doublet using subtomogram averaging from cryoelectron 
tomograms. We found no difference between the central and 
distal regions in neither the wild-type nor the mutant axonemes 
analyzed in this study, except for the presence or absence of the 
beaks. Therefore, what we refer to as the distal region in the 
following text includes the central region. We used the minus 
sign (— ) and plus sign (+) after the MTD number to indicate 
the proximal (which corresponds to the minus end) and cen- 
tral/distal (the plus end) regions, respectively. For example, 
MTD3+ means the central/distal part of MTD number 3. An 
MTD number without a sign means that the observed feature is 
uniform along the particular MTD. 

We first attempted to localize dyneins in the 96-nm unit in 
the central/distal region, which span almost 80% of the flagel- 
lum. Dyneins c, e, and / were localized in our previous study 
(Bui et al., 2008). However, there have been ambiguities be- 
tween dyneins a and d and between dyneins b and g because no 
mutants have been available that lack one of the pairs. Using 
idalO flagella, which have significantly reduced amounts of 
dyneins b, c, d, and e (Yamamoto et al., 2010), we were able to 




RSI RS2 

Figure 2. Dynein arrangement within the 96-nm unit in the distal region. 

(A] Comparison of the MTD2-8+ structures in wild type (WT) and idalO 
(which has significantly reduced amounts of dyneins b, c, d, and e). In 
wild-type MTD2-8+, a dynein density labeled with two names (such as 
a/d) is the one whose identity has been previously unsettled between 
the two species. The black arrows in idalO MTD2-8+ indicate the miss- 
ing densities that were determined to be those of dyneins b and d in this 
study. Bar, 1 6 nm. (B) Model of the common IDA architecture in the 96-nm 
unit showing the complete assignment of all major dynein isoforms. The 
coloring scheme is as follows: red, IDA; light blue, ODA; turquoise, RS; 
green, DRC; yellow, dynein f IC/LC. The same coloring scheme is used 
throughout all figures and supplemental figures unless mentioned other- 
wise. Surface renderings are always positioned with the proximal part on 
the left and the distal part on the right unless mentioned otherwise. 

assign dyneins a, b, d, and g in the averaged map of the central/ 
distal region (Fig. 2 A). This has completed the assignment of 
the common IDA architecture (Fig. 2 B). 

There are variations in dynein structures 
among MTDs and along the flagellum 

Based on our division of flagellar regions, we compared the 
dynein arrangements in different MTDs and in different re- 
gions. The averaged map of MTD2-9 indicates that the 96-nm 
unit structure is mostly preserved in the PEA region except for 
the absence of dyneins / and b and a significant reduction of 
ODA (Fig. 3 A). In this region, adjacent MTDs are linked by 
the circumferential interdoublet linkers (Fig. 1 C, red arrows; 
Ringo, 1967; Hoops and Witman, 1983). As a result of the low 
number of 96-nm repeating units present (two to three per 
MTD), dynein isoforms and the radial asymmetry were not 
determined in this region. 

Next, we focused on the structural differences among 
MTDs in the proximal and distal regions (Videos 2-4). We 
previously observed major differences in dynein distribution 
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Figure 3. Dynein structure and arrangement in 
the whole axoneme. (A) Dynein arrangement in the 
PEA region. Surface rendering of averaged map of 
MTD2-9 shows that dyneins b and f (black arrows) 
are missing. WT, wild type. (B) Comparison between 
the proximal and distal regions in the MTD groups 1, 
2-8, and 9. Minus (— ) and plus (+) signs denote the 
proximal and central/distal regions, respectively. The 
black arrows indicate the missing dynein b in MTD2- 
8—, MTD9— , and MTD9+. The orange arrowhead 
indicates the position of the 1-2 bridge between 
MTD1— and MTD2 — . The blue arrow points to the 
bulb density that differs between MTD1 — and MTD1 + 
(referred to as MD2 in Fig. 4). Structures shown in 
orange are the 1-2 bridges, and those in dark green 
are IDL3, an interdoublet linker between MTD1 and 2 
(Bui et al., 2009]. (C) View of the section cut perpen- 
dicular to the 1-2 bridge of MTD1 + and MTD1 -. The 
black arrowhead shows the dynein f IC/LC in MTD1 +, 
which is missing in MTD1 — . (D) Close-up views of a 
region proximal to RSI (black frame in B] of MTD1 — 
(left), MTD1+ (middle), and the overlapped image 
(right). The blue mesh in the right image represents the 
IDA density of MTD1 + corresponding to the region in 
red in the middle image. The blue arrowheads point to 
the tails of MD2 and dynein a. The white arrowheads 
indicate the bulb density on the A-tubule that is in con- 
tact with dynein fa and MD2. Bars, 16 nm. 
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among MTD groups 1, 2-8, and 9 in the distal region (Bui et al., 
2009). The dynein structures in MTD2-8-, MTD9 — , and 
MTD9+ are similar, all lacking dynein b (Fig. 3 B, black 
arrows; we have already reported that MTD9+ lacks dynein b; 
Bui et al., 2009). MTD1 — has a structure distinctively different 
from other MTDs and also from MTD1+ (Fig. 3 B). The most 
prominent difference is the presence of the 1-2 bridge (Fig. 3 B, 
orange arrowhead), as reported by Hoops and Witman (1983). 
MTD1— lacks dynein c (as in MTD1+; reported in Bui et al. 
[2009]) and the DHC and IC/LC of dynein / (blue arrow in 
Fig. 3 B). The 1-2 bridge and the IC/LC of dynein/would have 
a spatial conflict if they existed together (Fig. 3 C). The bulb at 
the position of dynein a in MTD1 — has a tail similar to that of 



dynein a (blue arrowheads in Fig. 3 D). However, the head of 
this bulb is bigger than that of dynein a and somehow elongated 
like dynein fa. This head is shifted proximally compared with 
the dynein a head (Fig. 3 D, right image). We interpret this 
structure as an MD candidate (MD2), as it has a ring-shaped 
head, characteristic of a DHC (discussed later; Fig. S2). 

The proximal region contains a distinct: set: 
of dyneins 

Having observed the differentiation in dynein arrangement, 
we wanted to pinpoint the position of respective dyneins 
in the proximal region, including those of MDs. To probe 
the presence of MDs in various dynein-deficient mutants, we 
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Figure 4. Localization of MDs. (A) Western blot 
analysis of wild-type (wt), ida4, ida5, idalO, and pf3 
axonemes using antibodies to DHC3, DHC4, DHC1 1, 
and IC140 (an IC of dynein f, used as a positive con- 
trol). The molecular masses of DHC3, DHC4, and 
DHC1 1 are expected values, as they migrate signifi- 
cantly more slowly than the 250-kD molecular mass 
marker, which is not shown. (B) The 96-nm surface 
renderings of MTD2-8+ and MTD1 + of wild type and 
MTD2-9- and MTD1 - of wild type, ida4, ida5, 
idalO, idal, and pf3. There are three MD densities 
that deviated from the common inner-arm dynein struc- 
tures, MD1-3, indicated with arrowheads in differ- 
ent colors. They can be the candidates of the MDs. 
A filled arrowhead indicates the presence of dynein, 
whereas an open arrowhead with the same border 
color indicates the absence of dynein from the respec- 
tive position. Bar, 16 nm. 



performed Western blot analysis on the ida4, ida5, idalO, 
and pf3 axonemes (analyses of the ida4 and ida5 axonemes 
have been already reported by Yagi et al. [2009]). The results 
showed that in the idalO axoneme, DHC3 and DHC4 are 
almost completely missing, whereas DHC11 is reduced 
but still detectable (Fig. 4 A). The pf3 axoneme shows a 
slight reduction in DHC4 and DHC1 1 and a significant reduc- 
tion in DHC3. 

Based on the Western blot results (Fig. 4 A and Table 1), 
we compared the structure of the flagellar proximal region of 
wild type, ida4 (lacking dyneins a, c, d, and DHC 11), ida5 (lack- 
ing dyneins a, c, d, e, and DHC 11 and most of DHC3 and 
DHC4), idalO (lacking DHC3 and DHC4 and partially lacking 



DHC1 1 and the DHCs of dynein b, c, d, and e), and pf3 (lacking 
various DRC subunits and dynein e and partially lacking DHC3 
and DHC 1 1 ; Figs. 4 B and 5). In addition, we analyzed the struc- 
ture of the idal axoneme, which lacks dynein/. 

There are several bulb densities that have shapes different 
from those of major species of dyneins in wild-type axonemes. 
These exceptional densities could be those of MDs. The first candi- 
date, referred to as MD1, is the bulb density between RSI and RS2 
in MTD1, which we previously suggested to be an MD (Bui et al., 
2009). The shape of this density appears similar between MTD1 — 
and MTD1+ (black arrowheads in Figs. 4 B and 5). Because its tail 
originates from the same position as the tail of dynein e (Fig. 4 B, 
wild-type image), MD1 could be an MD replacing dynein e in 
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Table 1 . Characteristics of mutants in this study 



Strain 
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UnL l l 
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uni /urn 

N\U \ , N\UZ 
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a, c, d, e b 


DHCl 1° 


DHC3 and DHC4 


MD2 


0/5/5/1 
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pf3 


e, DRC (partial)' 




DHC3 significantly reduced, 
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1 /6/4/0 








DHC4and DHC11 reduced 






Ida 10 


b, c, d, e reduced d 


DHC3, DHC4 


DHCl 1 reduced 


MD2, MD3 


0/12/6/2 


idal 


f 




DHC3, DHC4, DHCl 1 


MD1-3 


1/10/8/0 



"Determined in Yagi et al. (2009) and this study. 

b Kamiya etal., 1991; Kagami and Kamiya, 1992; Yagi et al., 2009. 

'Gardner et al., 1 994; Piperno et al., 1994. 

d Yamamoto et al., 2010. 

"Kamiya etal., 1 991 ; Myster et al., 1997. 



MTD1. MD1 is significantly reduced or missing in ida5 and idalO 
(open black arrowheads in Figs. 4 B and 5). 

The second candidate, MD2, is the density at the dynein a 
position in MTD1— (green arrowheads in Figs. 4 B and 5). 
MD2 exists in all the mutants analyzed in this study, including 
ida4 and ida5, which lack dynein a. Although MD2 looks simi- 
lar to the DHC of / the dynein /-deficient axoneme of idal 
retains MD2 (Figs. 4 B and 5). 

The third candidate, MD3, is a density that appears at the dy- 
nein d position in the proximal region (blue arrowheads in Figs. 4 B 
and 5). In the proximal region of the idalO axoneme, a bulge is 
present at the dynein d position, as in the wild-type axoneme (blue 
arrowheads in Figs. 4 B and 5), despite the fact that dynein d is 
significantly reduced in idalO, as estimated by SDS-PAGE and 
Western blot analyses (Fig. 2 A; Yamamoto et al., 2010). Thus, 
this proximal density may represent some MD replacing dynein 
d. However, it is also possible that it originates from the small 
amount of dynein d remaining in the idalO axoneme. 

In addition to these candidates for MD species, densities 
assigned to major dyneins may well be the sites for some MDs. 
Immunological studies using a newly raised antibody for dynein e 
DHC revealed that the content of this DHC in wild-type axonemes 
is ^50% lower in the proximal region than in the central and distal 
regions; this was confirmed by both immunofluorescence micros- 
copy (Fig. 6) and immunoblotting analyses of short and long fla- 
gella (Fig. S3 B). However, in our cryo-ET images, the structure of 
the e position in the proximal region is similar to that in the distal 
region. This observation suggests that some of the dynein e in the 
proximal region of wild-type axonemes is replaced by MDs. 

Each one of the aforementioned MD candidates has a 
ring-shaped head typical of dynein (Fig. S2). Although avail- 
able evidence has not allowed us to clearly determine specific 
loci of individual MDs, it is quite likely that MDs replace major 
dyneins in the proximal region and in some MTDs (see Discussion 
for more details; Fig. 7). The interpretation of MD candidates is 
summarized in Tables 2 and 3. 

The connections between ODA and IDA 
are asymmetrically distributed 

Our analysis showed no asymmetry in the ODA structure itself 
but revealed an asymmetry in the outer-inner dynein link (OIDL). 



For a better description of the OIDL, we numbered the ODA- 1-4 
from proximal to distal direction in the 96-nm unit, with ODA- 1 
connecting to dynein/IC/LC (Fig. 8 A). We found five OIDLs in 
total. OIDL1 (Unking the dynein/IC/LC and the tail of ODA-1) 
and OIDL3a (linking the DRC and the tail of ODA-3), which 
have been previously reported (Nicastro et al., 2006; Bui et al., 
2008), exist in all 96-nm repeating units, except where dynein/ 
is missing (MTD1— ). The three other newly found OIDLs are 
(1) OIDL2 connecting the distal part of the dynein/IC/LC to 
the tail and the 7 DHC of ODA-2, (2) OIDL3b projecting from 
the proximal part of the DRC toward the y DHC of ODA-3, and 
(3) OIDL4 bridging a protrusion distal to the DRC and the tail of 
ODA-4 (Figs. 8 A and S4). OIDL2 and OIDL3b are present only 
in MTD6-8- (Fig. 8 B), whereas OIDL4 is very thin and found 
only in MTD2- and MTD7- (Fig. S4). 

At a higher magnification, OIDL2 is a density protrud- 
ing from the distal end of dynein/IC/LC, whereas OIDL3b 
is a bulb density protruding directly from the DRC toward 
the 7 DHC of ODA-3 (Fig. 8 B, purple). Thus, OIDL2 can 
be an additional component of dynein/IC/LC in MTD6-8 — . 
Likewise, OIDL3b can be an additional component of DRC 
in MTD6-8-. 

The observed OIDL pattern roughly divides MTDs into 
two groups. One group (MTD6-8 — ) has at least three OIDLs in 
the proximal region, whereas the other group (MTD2-5— and 9—) 
has at least two OIDLs (Table 3). 

Discussion 

Chlamydamanms axoneme has a complex 
dynein arrangement 

We have completed the assignment of seven major dynein 
isoforms in the common IDA architecture (MTD2-8+) of the 
Chlamydomonas axoneme. As stated in the introduction, our 
previous studies had left an ambiguity regarding the distinction 
in loci between dyneins a and d as well as between dyneins b 
and g. The use of the idalO axoneme, which has significantly 
reduced amounts of dyneins b, c, d, and e, enabled us to unam- 
biguously determine their respective loci. Gardner et al. (1994) 
showed by biochemical analyses a decrease in the dynein b con- 
tent in the pfZ axonemes, which lack certain DRC components. 
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A 




Figure 5. Sections showing the dyneins in 
wild-type and mutant axonemes. (A-C) Illus- 
tration of the cross-section plane (transparent 
dissecting plane; A) in B and C in relation to 
the isosurface structure. Cross-sections of den- 
sity maps of wild type (WT) and mutants in the 
proximal (B) and distal (C) regions in wild type, 
ida5, ida4, idalO, idal, and pf3 are shown. 
The color of each arrowhead is the same as in 
Fig. 4. Bars, 32 nm. 
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Their result might seem to indicate that dynein b is located close 
to the DRC; however, our Western blot analysis of the pf2 and 
pf3 axonemes did not detect any significant decrease in the 
amount of dynein b (Fig. S5). The reason for the discrepancies 
is unknown, but we speculate that the short length of pf2 and pf3 
flagella might have affected the analyses, as dynein b is absent 
in the proximal region, as shown in this study. Heuser et al. 
(2009) reported that the axoneme of pf3 lacks ID As in the posi- 
tion of dynein e and g (IAD4 and IAD5 in their terminology). In 
our reconstruction of the pf3 axoneme, dynein e is also missing, 
but dynein g is present (Fig. 5). 

ODA-IDA connections have been shown for the first time 
to exhibit asymmetry and polarity. In MTD6-8-, ODA-IDA 
association is tighter than in other MTDs because at least three 
of the four ODAs in the 96-nm unit are connected to dynein 
/ IC/LC or the DRC (Fig. 8 B). This indicates that the coor- 
dination between ODA and IDA activities could be regulated 



differently in a subset of MTDs. Interestingly, the heads of the 
7 DHCs in ODA-2 and ODA-3 are in direct contact with protru- 
sions from the dynein/IC/LC (OIDL2) and the DRC (OIDL3b; 
Fig. 8 B). From these images, we are inclined to speculate that 
the ODA 7 DHC might function as a regulator that senses me- 
chanical strain. A mechanosensing role of the 7 DHC has also 
been proposed by Patel-King and King (2009) based on the ob- 
servation of the contact between the 7 DHC and the A-tubule. 
Similarly, dynein fa might also act as a mechanical strain sen- 
sor by contacting with a density blob on the A-tubule surface 
(Fig. 3 D, white arrowheads). 

The arrangement of the IDAs in MTD2-8+ is similar to 
that in the averaged structure of Tetrahymena cilia, sea urchin 
sperm flagella (Pigino et al., 2012), and mouse respiratory cilia 
(Ueno et al., 2012), with some minor deviations. Therefore, a 
common IDA architecture might be shared among species. On 
the other hand, the radial asymmetries of IDA in the nine MTDs 
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Figure 6. Immunofluorescence localization of inner-arm dynein e. 

(A) Immunoblot analysis of wild-type (wt), ida4, ida5, idalO, and pf3 
axonemes using antibodies against dynein e heavy chains and IC140 
(positive control). The molecular mass of dynein e is the expected value, 
as dynein e migrates significantly more slowly than the 250-kD molecular 
mass marker, which is not shown. The reduction of dynein e in idalO is 
consistent with the previous result (Yamamoto et al., 2010). (B) Indirect 
immunofluorescence observation of inner-arm dynein e (left; green signal) 
and a-tubulin (middle; purple signal) in wild-type and ida5 cells. The right 
images show the merged images. Note that almost no dynein e signal 
was detected in ida5 axonemes. The dynein e signals in the proximal 
region of wild-type axonemes are significantly reduced compared with 
the distal regions. The boundaries are indicated with white arrowheads. 
Bar, 10 pm. (C) Magnified images of part of a wild-type cell and the den- 
sitometric evaluation of the signal intensities of dynein e and acetylated 
tubulin along the flagellum. The nucleoflagellar apparatus (see Materials 
and methods) was analyzed. The entire length of flagella is shown by 
double-headed arrows. Bar, 5 pm. 



in Chlamydomonas flagella (Bui et al., 2009; this study) and 
mouse respiratory cilia (Ueno et al., 2012) indicate that each 
species might tailor their IDAs in a distinct manner correspond- 
ing to their specific functions. 

MDs are possibly localized to specific loci 
in the axoneme 

In addition to determining the arrangement of major dyneins, 
we have proposed the candidate loci for MDs (Table 2), 
although we have been unable to determine the specific rela- 
tionship between these loci and particular DHC species. 

In our Western blot, DHC3 and DHC4 are equally signifi- 
cantly reduced in the ida5 and idalO axonemes (Fig. 4 A). 



Therefore, it is not possible to distinguish between the DHC3 
and DHC4 loci from the results of Western blotting and cryo-ET. 
Western blot analyses of short and long flagella have suggested 
that both DHC3 and DHC4 localize proximally in wild-type 
axonemes (Fig. S3 ; Yagi et al., 2009). On the other hand, immuno- 
fluorescence microscope observations indicated that the amount 
of dynein e in wild-type axonemes was significantly reduced in 
the proximal region (Fig. 6), whereas our cryo-ET density maps 
did not show any reduction in the dynein e density in the 
proximal region of the wild-type, idal, or ida4 axonemes 
(Figs. 4 B and 5). Therefore, DHC3, DHC4, or both of them 
may well be replacing dynein e in the proximal region of these 
axonemes (Fig. 7). 

The relationship between MD1 and minor DHCs as well 
as between MD2 and minor DHCs is unclear. MD1 occurs all 
along the length of the MTD1. As only one species of MD, 
DHC 12, is known to localize along the entire axoneme length 
(unpublished data), this DHC might be MD1. However, our 
preliminary observations indicated that DHC 12, but not MD1, 
is greatly diminished in ida4. Therefore, the identity of MD1 as 
well as the localization of DHC 12 must await further studies. 

MD2 occurs at the dynein a site on MTD1 — and is present 
even in mutants that lack dynein a as a result of the loss of its 
subunit p28 or actin (ida4 and ida5). It clearly exists in ida5, 
which has significantly reduced amounts of DHC3, DHC4, and 
DHC1 1. However, as the amount of the dynein expected to local- 
ize to the MD2 site should be very small (^4/50 of major inner- 
arm dyneins), minor inner-arm dynein species that are detected 
only faintly in Western blot in ida.4 or ida.5 can be MD2. 

MD3 is the density that localizes at the dynein d position 
in the proximal region of all MTDs. It is absent in ida4 and ida.5 
and reduced in idalO (Fig. 3 B). On the other hand, DHC1 1 has 
been shown to be present only in the proximal region of the 
wild-type flagellum. It is absent in ida4 and ida5 and reduced in 
idalO (Fig. 4 A; Yagi et al., 2009). As DHC11 and MD3 show 
the same distribution pattern among mutants, DHC 11 might 
well be MD3. However, an analysis of dynein composition in 
short and long axonemes using ion-exchange chromatography 
showed that the amount of dynein d did not decrease in short 
flagella (unpublished data). Also, p38 and p44, the LCs of 
dynein d, have been shown to localize all along the length of the 
axoneme (Yamamoto et al., 2006, 2008). These observations 
suggest that dynein d is present from base to tip of the axoneme. 
Therefore, the postulated dynein d replacement by DHC11, if 
any, must occur only partially. We cannot exclude the possibil- 
ity that part or all of the DHC 1 1 molecules are replacing some 
other major dyneins in the proximal region of axoneme. 

In addition, there are a few peculiar dynein distribution 
patterns that we still cannot explain. In the axoneme of i da5, the 
bulb density at the dynein g position is significantly reduced in 
the MTD1 (Figs. 4 B and 5) and MTD5 (not depicted) along the 
whole length. This mutant expresses novel actin-like protein 
(Kato-Minoura et al., 1998) instead of actin. Inner-arm dynein g, 
whose actin subunit is replaced by novel actin-like protein, may 
be unable to bind to these two MTDs. 

It is somewhat puzzling that the proximal region of 
MTD2-8 in the axoneme of idalO lacks dynein a (Fig. 4 B). 
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Figure 7. Possible assignment of dynein locations. 

(A) Assignment of the potential location of the three MDs 
(DHC3, DHC4, and DHCl 1) in the proximal region. The 
arrowheads indicate possible locations: possibly, DHC3 
and DHC4 are at the yellow arrowhead, and DHCl 1 is at 
the gray arrowheads. Densities at the other arrowheads 
are not assigned. DHC3 may be present also on MTDl +. 
Bar, 1 6 nm. (B) The proposed distribution of all the major 
dyneins along the flagella. The colors of the bar indicating 
possible localization corresponds to the MD or unknown 
density pointed by the same arrowhead colors in A. The 
short bar with the minus sign indicates the proximal re- 
gion of the MTD, whereas the long bar with the plus sign 
indicates the central and distal region. 



This is because the reduction of dynein a in the idalO axoneme 
is only slight, although idalO is deficient in cytoplasmic preas- 
sembly of dynein (Yamamoto et al., 2010). It may be that this 
dynein is produced in a slightly reduced amount or in a slightly 
different structure so that it cannot bind to the proximal portion 
where dynein association sites are somewhat different from 
those in the distal region. 

Finally, an unexpected finding was that dynein b is miss- 
ing in the proximal region (Yagi et al., 2009) without being 
replaced by another dynein. This is in contrast with dyneins d 
and e, which we suggest are replaced by MDs at least partially. 
This specialty of dynein b might suggest its specific function. 
The distribution of all dyneins with possible MD replacement is 
summarized in Fig. 7 B. 

The presence of many DHCs is not unique to Chlamydo- 
monas, as there are many dynein genes found by genome analysis 



in other species such as in Tetrahymena (25 axonemal dyneins; 
Wilkes et al., 2008) and sea urchin sperm (15 axonemal 
dyneins; Morris et al., 2006). Thus, we might expect heteroge- 
neity in dynein composition along the length of cilia and fla- 
gella in other species as well. In human respiratory cilia, ODA 
is shown to be different in composition between the proximal 
and distal regions. DNAH5, an orthologue of the Chlamydomo- 
nas ODA y heavy chain, is present along the whole cilium, 
whereas DNAH9, an orthologue of the Chlamydomonas ODA 
(3 heavy chain, is found only in the distal two thirds of the cilium 
(Fliegauf etal.,2005). 

Proximal region of axoneme is highly 
differentiated from the main region 

As is clear from the differentiation of dynein species, the proximal 
region of the axoneme displays more extensive structural 
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Table 3. Summary of main features observed in the proximal and distal regions of Chlamydomonas axoneme in this study 
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Inner-arm dyneins in the distal region of the axoneme are arranged from proximal to distal location within the 96-nm repeating unit. Underscores indicate empty loci. 
Question marks indicate that dynein e might be replaced by an MD in the proximal region. 



heterogeneities among MTDs. This implies that the proximal 
region plays an important role distinct from that of the other 
axonemal regions: for example, a role for the bending initiation 
and waveform formation. The distal region of the flagellum 
might just follow the lead of the proximal region. 

The circumferential differentiation of the OIDLs divides 
the MTDs into several groups: two sides that undergo exten- 
sive shearing (i.e., MTD2-4 and MTD6-8) and those that do 
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Figure 8. Asymmetry and polar differences in the number of the OID 
linkers. (A) Surface rendering of MTD6-8- showing all the OIDLs (OIDLl, 
OIDL2, OIDL3a and b, and OIDL4) found between ODAs and IDAs in 
the axoneme. ODAs are numbered 1-4 from proximal to distal direction. 
ODA-1 is the one with the tail contacting to the dynein f IC/LC. DP is 
a density called the distal protrusion. OIDL4 is visible in not all MTD6- 
8 — ; it is present in specific MTDs such as MTD2— and MTD7— (Fig. S4). 
(B) Highlights of the differences between OIDLs in MTD2-9+, MTD2-5,9- 
and MTD6-8 — . From the comparison, OIDL2 and OIDL3b (colored in 
purple) are extra density protruding from the dynein f IC/LC and the DRC, 
respectively. The density of OIDL2 and OIDL3b was determined from the 
difference map between MTD6-8— and MTD2-4— . The proximal end of 
the doublet is toward the lower part of the image. Bar, 16 nm. 



not undergo extensive shearing (i.e., MTD5 and MTD9-1). 
MTD1 has structures clearly distinct from other MTDs, 
whereas MTD5 and MTD9 have connections to the DRC 
slightly different from other MTD groups (unpublished data). 
Compared with MTD2-4, MTD6-8 displays greater structural 
differences between the proximal and distal regions in the 
number of OIDL (Fig. 8 B). Therefore, there should be a func- 
tional difference between these two groups of MTDs (Fig. 9). 
In an organism like Chlamydomonas displaying asymmetric 
flagellar waveforms, sliding movements in MTD2-4 produce 
the force for the power stroke (or the principal bend; the bend 
that produces the largest force), whereas sliding in MTD6-8 is 
responsible for the recovery stroke (reverse bend), which 
requires less force than the power stroke. The asymmetry in 
the OIDLs among MTDs may be related to such asymmetry in 
waveform or force production properties. 

Finally, there must be a mechanism that determines the 
complex arrangement of dyneins in the axoneme. As we ob- 
served longitudinal differences in the arrangement of micro- 
tubule inside protein la (MlPla; Pigino et al., 2012), which is in 
the vicinity of dynein tails (Fig. 4; Pigino et al., 2012), we spec- 
ulate that MlPla provides the platform for the assembly of IDA 
and other flagellar components and functions as the molecular 
ruler to partition the proximal, central, and distal regions. 




Figure 9. Summary of the observed structural polarization and asymme- 
try in the axoneme and their possible importance for bending motion. The 

division of the axoneme into MTD groups based on the OIDLs structures. 
The group of MTD2-3-4 (dark green] has almost the same OIDLs in the 
proximal and central regions, whereas the group of MTD6-7-8 (blue) has 
extra OIDLs and density in the proximal region than in central regions. 
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Conclusion 

We showed in this study intricate radial and longitudinal vari- 
ations in the Chlamydomonas axoneme. Which feature is uni- 
versally conserved, and which is specific to Chlamydomonas! 
These questions can be answered by systematic structural 
studies of axonemes in other organisms. We already know 
that some of the features observed, such as the 1-2 bridge, 
are specific to Chlamydomonas. On the other hand, other fla- 
gella have their own unique features such as the 5-6 bridge 
in the flagella of sea urchin sperm or the parafiagellar rod in 
Trypanosome. It seems clear that the eukaryotic flagellum has 
a highly elaborate structure that can be differentiated longitu- 
dinally and circumferentially. 

Our present structural study poses an exciting challenge 
for future studies of flagellar function. New assays must be 
developed to explore the function of a component present only 
in specific regions of a subset of MTDs. Elucidation of the 
biophysical properties of MDs will provide important clues to 
the functional significance of these dyneins and the axonemal 
proximal regions. 

Materials and methods 

Cell culture and flagella purification 

Chlamydomonas reinhardtii wild-type 137c and mutants were cultured 
using tris-acetate-phosphate medium (Gorman and tevine, 1965). 
Deflagellation was induced by the dibucaine method of Witman 
(1986]. The separated flagella were centrifuged at 5,200 g for 40 min 
at 4°C. The pellet was suspended in HMDEKP solution (30 mM Hepes, 
5 mM MgS0 4 , 1 mM DTT, 0.5 mM EDTA, 25 mM KCI, and 0.5% 
polyethylene glycol [molecular weight of 20,000], pH 7.4) at ice-cold 
temperature. The flagella are demembranated by HMDEKP + 0.5% NP- 
40 (EMD Millipore), spun down again at 5,200 g for 40 min, and 
resuspended in HMDEKP. 

SDS-PAGE and Western blot 

The composition of DHCs was analyzed by SDS-PAGE with a 4% poly- 
acrylamide and 6 M urea gel. 10 pg of axoneme was applied on the gel 
for each strain. Immunoblot was performed by the method of King et al. 
(1986). pAbs produced against the N-terminal polypeptides of DHC3, 
DHC4, DHC5 (dynein b), and DHC1 1 were used (Yagi et al., 2009). In 
addition, a pAb against the N-terminal sequence of DHC8 (dynein e) cor- 
responding to amino acids 21-714 was produced in this study. Antibody 
against IC 1 40, an IC of inner-arm dynein f, was also used (a gift from 
W. Sale, Emory University, Atlanta, GA). Immunoreactive bands were 
detected using HRP-conjugated secondary antibody (GE Healthcare) and 
chemiluminescence reagents Chemi-Lumi One Super (Nacalai Tesque, Inc.) 
or ImmunoStar LD (WAKO Chemicals USA). 

Immunofluorescence microscopy 

Immunofluorescence microscopy was performed according to Taillon 
etal. (1992). Wild-type and ida5 cells were demembranated in 1% NP-40, 
fixed with 2% PFA for 2 min at room temperature, and treated with — 20°C 
acetone. Nucleoflagellar apparatuses (complex of two axonemes, two 
basal bodies, and a nucleus) produced by the demembranation of CW92 
cells were also used (Wright et al., 1 985). The fixed samples were stained 
with an a-tubulin mAb (T5168; Sigma-Aldrich) and with a specific pAb of 
dynein e. Alexa Fluor 488-labeled anti-rabbit IgG (Invitrogen) or rhodamine- 
labeled anti-mouse IgG (Sigma-Aldrich) was used as the secondary anti- 
body. Images were acquired using an 1X70 microscope with a 1 00x/l .3 
NA objective lens (Olympus) and a cooled charge-coupled device camera 
(ORCA-AR; Hamamatsu Photonics). 

EM 

The demembranated flagella were quick frozen on holey carbon grids 
(R2/1 ; QUANTIFOIt) using the Vitrobot plunger machine (FEI). 3.5 pi of 
10-nm gold colloidal solution (Sigma-Aldrich) was deposited twice on the 



grid to use as the fiducial markers in tomogram reconstruction. We col- 
lected tilt series in the range of ±60 degrees with a 2-degree increment 
at 19,303x magnification at cryotemperature (approximately — 175°C) 
using a transmission electron microscope (200 kV; Tecnai F20; FEI) 
equipped with a Tridiem energy filter (operated at 20eV slit; Gatan, Inc.) 
and an UltraScan 1000 camera (2,048 x 2,048 pixels; Gatan, Inc.) 
with defocus ranging from 4 to 6 urn. Data of idal mutants were col- 
lected using a transmission electron microscope (300 kV; Tecnai Polara; 
FEI) equipped with a GIF2002 energy filter (Gatan, Inc.) and a MultiScan 
camera (2,048 x 2,048 pixels; Gatan, Inc.) with defocus ranging from 
8 to 9 pm. In both microscopes, the tilt series acquisition was controlled 
using Explore 3D (FEI). The pixel size of the tilt series is 7.25 A. 

The size of our camera prevents us from analyzing the whole flagel- 
lum simultaneously. With our acquisition setup, a tomogram contains only 
a region of ~1 .4 pm, whereas the length of the whole flagellum is normally 
between 10 and 15 pm. Therefore, we systematically collected the tomo- 
grams of different regions along the flagellum by acquiring the tilt series 
and the corresponding low-magnification image for later identification. To 
increase the chance of proximal region of the axoneme during acquisition, 
we identified the proximal region by its square-end microtubule compared 
with the hanging-end microtubule of the tip region. We paid attention to 
collect noncompressed flagella by acquiring tilt series of axonemes that 
had diameter of ~220 nm. 

Subtomogram averaging 

After the tomographic reconstruction by IMOD (Mastronarde, 1997), the 
flagella were classified into different regions by visually examining their 
cross-sections according to the region-specific features such as the 1-2 
bridge (proximal), beak (central), or no beak (distal), as shown in Fig. 1 B. 
As the PEA region contained few repeat units and coexisted with the proxi- 
mal region, subtomograms belonging to the PEA regions were aligned 
together with the proximal region and separated before averaging. 

Subtomograms of 200 x 200 x 200 pixel (~1 44 nm in each dimen- 
sion) were extracted with roughly 96-nm periodicity from the same region 
of the flagella by Bsoft (Heymann, 2001 ) and then gone through the 96-nm 
repeat alignment separately, as described elsewhere (Ishikawa et al., 
2007; Bui et al., 2008, 2009), using a combination of software pack- 
ages, including SPIDER (Frank etal., 1996), Bsoft (Heymann, 2001), TOM 
(Nickell et al., 2005), and an in-house program. In short, the subtomo- 
grams from the same MTD from one tomogram were aligned along the mi- 
crotubule based on cross-correlation (intradoublet alignment) to make 
subaverages along each MTD. The obtained subaverages were then 
aligned between MTDs. 

After alignment, the position of the doublet in the entire axoneme 
was determined using the specific features of Chlamydomonas (Fig. 1). 
Subtomograms belonging to each MTD and region were then averaged 
together, refined with the new average and corrected for the missing 
wedge to ensure isotropicity. Isosurface visualization was performed in 
UCSF Chimera (Pettersen etal., 2004). The segmentation (coloring) of the 
maps for outer- and inner-arm dyneins, RSs, the DRC, and dynein f IC/tC 
was performed using the UCSF Chimera program and based on previous 
studies by ourselves (Ishikawa et al., 2007; Bui et al., 2008; Pigino et al., 
201 1, 2012) and others (Nicastro etal., 2006; Heuser et al., 2009). The 
96-nm model in Fig. 2 B was created by tracing density map by 3DMOD 
(Kremer et al., 1996). 

To quantify the distortion by compression, we fitted an ellipse on the 
model of the middle cross-section of the flagella built from the alignment 
parameters and defined the elliptical distortion as the long axis divided by 
the short axis of the fitted ellipse. Only flagella with elliptical distortion 
<20% were included in the average (Fig. SI A). By Fourier shell correla- 
tion, the resolution of the averaged map of MTD2-8— at 0.5 cut-off was 
36.8 A (Fig. SI B). 

To ensure that the structural features that we observed were repro- 
ducible, we split our wild-type data into two sets and examined whether 
the averages from those two sets showed similar structural features. In 
addition, we also averaged the mutants used in this study to double check 
the features that were not missing in the mutant were similar to those of 
wild type. 

Density maps 

The density maps reconstructed in this study were deposited in EMDataBank 
as EMD-21 13 to 21 18 and EMD-2121 to 2123 (MTD 1-9 of the distal re- 
gion from wild type), EMD-21 19, 2020, and 2 124 to 2130 (MTD 1-9 of the 
proximal region from wild type), and EMD-21 32 and 2131 (average of the 
distal and the proximal regions from wild type, respectively). 
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Online supplemental material 

Fig. SI describes the quality of our image analysis. Fig. S2 shows struc- 
ture of dynein heads (AAA rings) as sections from averaged tomograms. 
Fig. S3 shows immunoblot analysis of inner-arm dyneins in long and short 
flagella. Fig. S4 demonstrates connections between IDAs and ODAs. 
Fig. S5 shows Western blot of dynein b in wild-type, pf2, pf3, and idalO 
axonemes. Video 1 shows structure of axoneme at the proximal, central, 
and distal regions. Video 2 highlights the difference between the proxi- 
mal and the central/distal regions in MTD2-8. Video 3 demonstrates the 
difference between MTD2-4 and MTD6-8 in the proximal region. Video 4 
shows the difference between the proximal and the distal regions in 
MTD1 . Online supplemental material is available at http://www.jcb.org/ 
cgi/content/full/jcb.201 201 1 20/DC1 . 
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